LINEARIZED PLASTICITY IS THE EVOLUTIONARY 
r-LIMIT OF FINITE PLASTICITY 



ALEXANDER MIELKE AND ULISSE STEFANELLI 



Abstract. We provide a rigorous justification of the classical linearization 
approach in plasticity. By taking the small-deformations limit, we prove via 
F-convergence for rate-independent processes that energetic solutions of the 
quasi-static finite-strain elastoplasticity system converge to the unique strong 
solution of linearized elastoplasticity. 



1. Introduction 

This paper is devoted to the rigorous justification of the classical linearization 
approach in finite-strain elastoplasticity. When restricting to the small-deformation 
realm it is indeed customary to leave the nonlinear finite-strain frame and resort 
to linearized theories instead. This reduction is usually motivated by means of 
heuristic Taylor expansion arguments. Here, we aim at complement these formal 
motivations by providing a rigorous linearization proof by means of an evolutionary 
T-convergence analysis of rate-independent processes. In particular, we address the 
general time-dependent case, which e.g. allows for cyclic loading. 

In the stationary framework, the pioneering contribution in this context goes 
back to Dal Maso, Negri, & Percivale [DNP02] who devised a convergence 
proof of finite-strain elasticity to linearized elasticity. Later, the argument has 
been refined by Agostiniani, Dal Maso, & DeSimone [ADDll] and extended 
to multi-well energies by Schmidt [Sch08] and to residually stressed materials by 
Paroni & Tomassetti [PT09, PTll]. The reader is also referred to [GNIO, MNll, 
NeulO] for some related results in the direction of homogenization, to [ADll] for an 
application to the study of nematic elastomers, to [BSV07, Sch09] in the context of 
convergence of atomistic models, and to [SZll] in relation with dislocation theory. 

To our knowledge, this is the first result in the evolutionary case. With respect 
to the stationary case of [DNP02], the evolution situation is quite more involved. 
Indeed, the argument in [DNP02] relies on the F-convergence proof of the small- 
deformation energy functional to its linearization limit. Here, we are instead forced 
to cope with the occurrence of dissipative plastic evolution by means of a deli- 
cate recovery sequence construction relating energy and dissipation. We emphasize 
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that finite-strain elastoplasticity is based on the multiphcative decomposition of 
the strain tensors. Moreover, the plastic tensor is to be considered as an element 
of a multiplicative matrix group. We have to control these noncommutative multi- 
plicative structures in linear function spaces and to establish their convergence to 
the corresponding linear additive structures. In order to give some details in this 
direction wc cannot avoid introducing some minimal notation. 

Finite-strain elastoplasticity is usually based on the multiplicative decomposition 
Vip = FciFpi [Lee69]. Here : O R** is the deformation of the body with respect 
to the reference configuration C M.'^ {d = 2, 3) while Fe\ and Fpi € SL((i) stand 
for the elastic and the plastic strain, respectively. Then, the stored energy in the 
body is written as 



where Wei is a frame-indifferent elastic stored-energy density and Wh describes 
hardening. The plastic flow rule is expressed by means of a suitably defined dissi- 
pation distance D : SL{d) x SL{d) — >■ [0, oo]. In particular D{Fpi,Fpi) represents 
the minimal dissipated energy for an evolution from the plastic strain Fpi to Fpi 
and is given via a positively 1-homogcncous dissipation function R by 



the infimum being taken among all smooth trajectories P : [0, 1] M.'^^'^ connecting 
-Fpi to _Fpi. Staring from these functionals, by specifying loadings, boundary, and 
initial conditions, suitably weak solutions of the quasi-static finite-plasticity system 
(see Section 2) can be defined. We refer to [Mie03] for more information on the 
mathematical modeling of finite-strain elastoplasticity. There also models with 
additional hardening variables are given. Here we however refrain from maximal 
generality in order to emphasize the main features of the limiting process. 

Let now the deformation and the plastic strain be small. In particular, for e > 

let = id+e?i and Fpi,e = I+ez where u is interpreted as the displacement and z 
is the linearized plastic strain. Correspondingly, we have that i^ei,e = ^f£-^pi\ — 
(id-|-£Vu)(/-|-£2;)~^ and we are lead to the consideration of the rescaled finite-strain 
elastoplasticity energy and dissipation functionals 



Note that the rescalings above are such that, by assuming Wei and Wh to admit a 
quadratic expansion around identity, one can check that 






-D((/+£Z),(/+£Z)). 
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This pointwise convergence is the classical justification of linearization in plasticity. 
On the other hand, it is not sufficient in itself for proving that finite-strain elasto- 
plasticity trajectories actually convergence to a solution of the linearized-plasticity 
system. 

Before going on let us mention that the solution concept which is here under 
consideration is that of energetic solutions. Starting from [MT04], this solution 
notion has been extensively applied in many dififerent rate-independent contexts. 
We shall however record that one of the main motivations for introducing energetic 
solutions was exactly that of targeting existence theories for finite-strain clastoplas- 
ticity. In this respect, note that the only available existence result for finite-strain 
elastoplastic evolution has been recently obtained within the energetic solvability 
frame in [MM09] after adding the regularizing term |VFpi|'' for r > 1 (see also 
[MM06] for some preliminary result). 

Our result consists in proving the convergence of energetic solutions of the finite- 
strain elastoplasticity system to linearized-plasticity solutions. In order to prove 
this convergence we follow the abstract evolutionary F-convcrgcnce theory for ener- 
getic solutions of rate-independent processes developed in [MRS08] . We shall men- 
tion that this evolutionary F-convergence method has recently attracted attention 
and has been successfully considered in connection with numerical approximations 
[KMR05, MR09, GP06a], damage [BRM09, TMIO], fracture [GP06b], delamina- 
tion [RSZ09], dimension reduction [FPZIO, LMll], homogenization [Tim09], and 
optimal control [Rin08, Rin09]. 

According to [MRS08], the convergence of the trajectories {u^,z^) follows by 
proving two separate F-liminf inequalities for energy and dissipation and con- 
structing of a mutual recovery sequence relating both. Note that separate F- 
convcrgcncc for energy and dissipation is not sufficient to pass to the limit within 
rate-independent processes. Apart from the additional technicalities due to the 
presence of the plastic strain and the dissipation functional, it is the delicate con- 
struction of the mutual recovery sequence that distinguishes our argument from all 
the already developed stationary analyses in the spirit of [DNP02]. 

2. Problem setup and results 

Let the reference configuration c M'' be an open, bounded, and connected set 
with Lipschitz boundary. Moreover, let F c dQ be relatively open with 'H'^~^(F) > 
0. We define the state space as 

Q:=UxZ:= {u € R^niW^) | u = on F} x L'^{n-W^'"^). 

Note that the choice of the homogeneous Dirichlet condition on the displacement 
u is just motivated by the sake of simplicity. In particular, different boundary 

conditions may be considered as well. 

For all given A G M.'^^'^ we denote its symmetric and antisymmetric parts as 
^sym ._ (^A+A^)/2 and A™" = ^-A^^". We indicate by R^^^ and Mf^^f sub- 
spaces of symmetric and antisymmetric tensors, respectively, whereas M^^J^ stands 
for the subspace of symmetric and trace-free tensors, also called deviatoric tensors. 
The standard Euclidian tensor norm is denoted by | • | and, for all A e M''^'* and 
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T > 0, Br{A) indicates the ball Br{A) := {B 6 R'^^'' | |^ - fi| < r}. Moreover, 
the symbol | • |t stands for the seminorm 

\A\l := ^A:TA 

where the 4-tensor T e R'^x'^x'^x'' jg symmetric (Tjj/j^ = T^^jj) and positive semi- 
definite. For finite-strain elastoplasticity we use the classical notations 

SL(rf) := {P e R^""^ I detP = 1}, 

SO{d) := {R e SL{d) \ R^R = RR'^ = I}, 

GL+{d) := {Q G R'^'"^ I detg > 0}. 

We assume that the elastic energy density functional Wei fulfills 

Wei : M.'''"' [0,oo], Wei G C^(GL+(rf)), Wei = 00 on M.'^'"' \ GL+{d), (2.1a) 

VF G GL+(d) VP G SO{d) : Wei(PP) = Wei(P), (2.1b) 

VP G GL+(d) : Wei(P) > cidist2(P,S0(d)), (2.1c) 

VP G GL+(d) : |P^9^^Wel(P)| < C2(Wel(P) + C3), (2.1d) 

aC > V5 > 3cei(^) > V^ G Pe.,(5)(0) : |Wei(/+A) - l^l^l < S\A\l, (2.1e) 

for some positive Ci, C2. Assumption (2.1b) is nothing but frame indifference and 
the nondegeneracy requirement (2.1c) is quite classical. Assumption (2. Id) entails 
the controllability of the Mandel tensor dpWciiF) by means of the energy. This 
is a crucial condition in finite-strain elastoplasticity (cf. [Bal84b, Bal02]) and was 
used in the context of rate-independent processes in [FM06, MM09]. Condition, 
(2.1c) encodes the local quadratic character of Wei around identity. More precisely, 
(2.1e) states that | • |c is the second order Taylor expansion of Wei at I, and may 
be reformulated by saying that A i->- Wei(/+A) is locally restrained between two 
multiples of | • namely, 

V^ > VA G Pee.w(0) : {i-S)\A\l < Wei(7+^) < {l+6)\A\l 
Moreover, (2.1e) entails 

Wel(/)=0, afWel(/)=0, a2Wel(/)=C, (2.2) 

which, in particular, yields that the reference state is stress free. On the other 
hand, by assuming (2.2) and letting Wei G in neighborhood of /, relation (2.1e) 
follows. 

Note that the symmetry of the elastic tensor C (implicitly assumed in the no- 
tation I • |c) may be directly obtained from the last of (2.2) by assuming addi- 
tional smoothness on Wei. Moreover, letting A G R'^^'^ be given, as we have 
that exp(A™*') G SO{d), the frame indifference (2.1b) entails that the function 
t 9j?Wei(exp(fj4''"*')) is constantly equal to dpWeiil) = 0. Hence, by taking its 
derivative with respect to t and evaluating it at f = we get Cj4'^"*' = 0. Namely, 
C necessarily fulfills also the so called minor symmetries Ciju = ^jiu = Cijek and 
we have 

y^gj^dxd. CA = CA^y"'. (2.3) 

On the other hand, as effect of the nondegeneracy (2.1c) and assumption (2.1e) 
we have that C is positive definite on Mf^. Indeed, by linearizing d{-,SO{d)) 
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around identity we have [FJM02, (3.21)] 

y^^jgdxd. d{B,SO{d)) ^\B'y'^-~I\+0{\B-I\^). (2.4) 
Hence, given A G R''^'^ and rj, S > 0, hy choosing B = I+rjA in the latter we have 

^ ,^sym|2 2L(f(I+rjA,S0(d)) 
rj— T] 

(2.1c) 1 (2.1o) 

< \im ^W^i{I+r]A) < {l+d)\A\l 

77— >0 T] 

SO that, by taking 5 — >^ 0, we have 

g M'^xrf : ci\A^^"'f < \A\l = lA^^'^ll. (2.5) 

Note that all assumptions (2.1a)-(2.1e) are consistent with the usual polycon- 
vexity framework 

F !->■ Wei{F) polyconvex, 
Wei(-F)-^oo for detF-)-0. 

Our assumptions on the hardening functional : K''^'^ — > [0, 00] read 
•^^^^^-too if PeM<^>«^\i^, 

where K is compact in SL(d) and contains a neighborhood of I, (2.6b) 

Wh : M'*^'* ^ M is locally Lipschitz continuous and (2.6c) 

3e > > 3ch(<5) > e Bch(5)(0) : | Wh(/+A) - \A\l\ < 5\A\l, (2.6d) 

3c3 > VA G M'*^'^ : Wh(-f+A) > C3\Af. (2.6e) 
Note that by assumption (2.6b) we can find a constant ck > such that 

P&K ^ \P\ + \P-^\<CK, (2.7) 

P&Sh{d)\K \p-l\>—, (2.8) 

Ck 

The rather strong technical assumption on Wh that its effective domain K = {P € 
SL{d) I Wh{P) < 00} fulfills (2.7) is crucial as it will provide L°°-bounds that are 
essential in order to control the multiplicative terms {I+eVu){I+ez)~^ . Moreover, 
by combining (2.6d) and (2.6e) we check that 

y^gj^dxd. C3|^P<|^|^. (2.9) 

As for the dissipation we assume that 

,00) convex and positively 1-homogeneous, (2.10a) 
VPGMlj^xj^: C4|P| < i?'^" (P) < C5IPI, (2.10b) 

R : M'^x'^ ^ [0, 00]; R{z) := ( ^""^^^ ^ ^ ^^ev > 

[00 else. 



(2.10c) 
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for positive C4, C5. Moreover, we define 




yP,QeKc SL{d) : D{P,Q) < cq, D{I,P) < ce\P-I\. 



(2.12) 



For the first estimate the continuity of D and the compactness of K is sufficient. 
For the second, we need to establish the estimate only for P close to /, where 
it follows from D{I,P) < R'^'''' {log P) < C5|logP| < ce\P-I\, since the matrix 
logarithm is well-defined and Lipschitz continuous in a neighborhood of I. See also 
[MM09, Ex. 3.2] and the references given there for global bounds on D. 

The quasistatic evolution of the finite-strain and linearized elastoplasticity sys- 
tems are driven by the energy functionals Wg, Wo : (—00, 00] given by 



Note that, if the second integral in the definition of >Ve(M, z) is finite, then I+ez G K 
almost everywhere by (2.6a). Hence, the inverse {I+ez)~^ exists and the first 
integral is well defined. 

We prescribe the generalized loading as 



and, by letting := e£, we introduce some notation for the total energy functionals 
£e, £0 : [0,T] X Q (-00, 00] as 

£e{t,U,z) ■.= We{u,z) --%{t),u) =We{u,z) - {m,u) 

£ 

£o{t,u,z) :=Wn{u,z)~ {£{t),u), 

Eventually, the dissipative character of the evolution is encoded into the dissipation 
functions Z>e, -Do : W^""^ xW^""^ [0,oo] and functionals P^, Do : {L^{n;R'^'"^)f 
[0, 00] given by 

Ds{zi,Z2) := -D{I+ezi,I+ez2), Do{zi,Z2) := R{z2-zi), 




eeW^''^{0,T;U') 



(2.13) 



e 



'Ds{zi,Z2) := / Ds{zi,Z2)dx, T>q{zi, Z2) := / £'0(21, 22) dx. 
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The total dissipation of the process over the time interval [0, t] C [0, T] will be given 

by 

f ^ 

Diss^,,(z; [0,t]) := sup <^ ^Pe(z(f ), z{f-')) \ {0 = t° < ■ ■ ■ < t"" = t} 

where the sup is taken over all partitions of [0,t\. 

From here on, we term Rate- Independent System (RIS) the triple {Q,£e,T^s) 

given by the choice of the state space Q and the energy and dissipation functionals 
£g and V^. The term evolutionary T-convergence refers to a suitable notion of con- 
vergence for rate-independent systems in the spirit of [MRS08] which in particular 
entails the convergence of the respective energetic solutions. 

A crucial structure in the energetic formulation of RIS is the set Sg{t) of stable 
states at time t £ [0,T], which is defined via 

Ss{t) := ^{u,z) G Q I £e{t,u,z) < oo and 

£eit,u,z) < £eit,u,z) +Ve{z,z) y{u,z) e q}. 

Our assumption on the initial data reads 

SM B {ul,z',) ^ K,zO) weakly in Q, zJJ e L'{n;RZ^), 

£,{0,Ue,Ze) ^So{0, u°,z'o)- (2-14) 

Note that the latter assumption is not empty as it is fulfilled at least by the natural 
choice {uo, zq) = (0, 0) if ^(0) = 0. 

Definition 2.1 (Energetic solutions). Let e > 0. We say that a trajectory : 
[0,T] — > {ue,Ze) G Q is an energetic solution (related to the RIS {Q,£e,'Dg)) if 
{ue{0), Ze{0)) = (m",z°), the map 1 1— ?■ {£,Ue) is integrable, and, for all t G [0,T], 

iu,{t),zS))<^S,{t), (2.15) 
£S,Mt),zS)) + 'DissvAze;[0,t])=£e{0,u°,z°)- [ {£,u,)ds. (2.16) 

^0 

An energetic solution will be called a finite-plasticity solution if e > and a 
linearized-plasticity solution for £ = 0. 

Note that linearized-plasticity solutions (wq, zq) are unique as effect of the qua- 
dratic and uniformly convex character of Wq. Moreover, from assumption (2.13) 
we get that (wo,^o) e W^'^{0,T; Q) and 

yt€[0,T]: Dissp,(2o;[0,i]) = / R{zo)ds. 

Jo 

The reader is referred to [Hil50, Lub90, Mar75] for some general introduction to 
plasticity and to [HR99, Joh76, SuqSl] for the classical well-posedness theory for 

linearized clastoplasticity. 

Our main result reads as follows and will be proved in Section 3 as a special 
instance of the general theory of [MRS08] . 
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Theorem 2.2 (Finite plasticity r-convorgos to linearized plasticity). Assume (2.1)- 
(2.6), (2.10), and (2.13)-(2.14). Let {us.Zs) he a finite-plasticity solution. Then, 
{us{t), Ze{t)) — >■ {uo{t) , zo{t)) weakly in Q for all t € [0,T] where (uo,^o) is the 
unique linearized-plasticity solution. 

Theorem 2.2 is exclusively a convergence result. In particular, we assume that 
finite- plasticity solutions exist. Note however that the existence of finite-plasticity 
solutions is presently not known within our minimal assumption frame. A pos- 
sibility here would be that of considering directly some more regular situations 
including extra compactifying terms like |VFpi|'' (r > 1) such that finite-plasticity 
solutions exist [MM09]. We shall not follow this line here but rather present a 
second result based on approximate minimizers of the related incremental prob- 
lems. Indeed, given the time partitions {0 = tl < ■ ■ ■ < t^' = T} with diameters 
Te := maxj=i^,..^/yf^(t^ — tl~^) — )• as £ — > 0, the (iterative) incremental problem 

{ui,zi)GAT:gmin{£,{ti,u,z) + V,{zi-\z)) for i = l,...,N, 

(«,t))GQ 

may not be solvable (cf. [CHM02], still see [Mie04, MM06] for some additional 

discussion). Hence, following [MRS08, Sec. 4] we fix a sequence < — ^ in 
order to control the tolerances for the minimizations and consider the following 
approximate incremental problem 

Find itcrativcly {v.l,zl) G Q such that 

£,(tl,vi,zl)+V,iz^\zl) (2.17) 
< {tl - tl~^)a, + inf(„,„)eQ {£e{tl,u,z) + V,{zl-\z)). 

By the definition of infimum the latter always admits solutions and we will show 
the following convergence result. 

Theorem 2.3 (Convergence of approximate incremental minimizers). Under the 
assumptions of Theorem 2.2 let (u*,^*) he approximate incremental minimizers 
and {ue,'Ze) he the corresponding right- continuous, piecewise- constant interpolants 
on the time partitions. Then, (u^{t),Jg{t)) {uo{t), zo{t)) weakly in Q for all 
t G [0,T] where (uo,zo) is the unique linearized-plasticity solution. 

In the finite-elasticity case (stationary), using ideas from [DNP02] the conver- 
gence of approximate minimizers has been considered in [PT09] . 

3. Proofs 

The argument basically follows the lines of the abstract analysis of [MRS08]. 
Still, our setting cannot be completely recovered from the application of the above- 
mentioned abstract theory as extra care is needed for the treatment of the multi- 
plicative nonlinearities. We hence resort in providing here an independent proof. 
After establishing the coercivity of the energy in Subsection 3.1, the proof strategy 
relies in providing two separate F-liminf inequalities for £e and Vg and a mutual 
recovery sequence argument relating both. This is done in Subsections 3.2 and 3.3 
below. Eventually, the proofs of Theorems 2.2 and 2.3 are outlined in Subsections 
3.4 and 3.5, respectively. 

A caveat on notation: henceforth the symbol c stands for any positive constant 
independent of s and 5 but possibly depending on the fixed data. In particular, 
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note that c may change from hne to hne. Moreover, in the foUowing we use the 
short-hand notation, for all A gR'^^'^, 

W^iiA) := lwei(/+eA), Wl{A) := lwh(J+£A), Wl{A) := ^W],{I+eA). 



3.1. Energy coercivity. We start by providing a uniform coercivity result for the 
energy. It follows the ideas in [DNP02] and relies on the Rigidity Lemma [FJM02, 
Thm. 3.1]. 

Lemma 3.1 (Coercivity). There exists c > such that, for all {u,z) e Q 

\\yu\\l. + \\Z\\1. + < C{l+We{u,z)). (3.1) 

Proof. Let us assume with no loss of generality that Weiu, z) < oo, so that I+ez € 
K almost everywhere by assumption (2.6a). Hence, |/+e2| < ck almost everywhere 
from property (2.7) and the inverse {I+ez)~^ exists almost everywhere. Thus, we 
have that ||£i;||L~ < c. Moreover, one readily checks from the coercivity (2.6e) that 

cs\\z\\l.< I W^{z)dx<We{u,z). (3.2) 
Jo, 

For the displacement u wo follow ideas from [DNP02]. Given any Q G SO(rf) by 
letting ^p = id+£u and fei = V<p(/+£2;)~^ we have 

|V^-Q|2 = |V(^ - Q{I+ez) + eQz\^ = \{F,i-Q){I+sz) + sQz\^ 
< c(|Fei-Qn/+e^|2 + e^\z\^) < c{\F,,-Q\^ + e^\z\^). 

In particular, by passing to the infimum for Q G SO{d) we have checked that 

dist^(V</j,SO(d)) < c(dist^(Fei,SO(rf)) + £2|^|2). 

By taking the integral in space and using the nondegeneracy condition (2.1c) we 
obtain that 

dist^(Vv>, SO(rf)) dx<c f dist^(Fei, SO {d)) dx + cs^ [ \z\^ dx 
Jn Jn 

(3.2) 

< s'^cWs{u,z). 
Hence, the Rigidity Lemma [FJM02, Thm. 3.1] ensures that 

\\Vip-Q\\l, <e\We{u,z) 

for some constant rotation Q € SO{d). Finally, using [DNP02, Prop. 3.4] and 
<^|r = id as u gU, we conclude |Q— < £^cWe(w, z). Then, we have 

W\h = ^l|V<^-/||^. < |||V^-g||£. + 1\\Q-I\\l, < cW,{u,z) 
and the bound (3.1) follows. □ 
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3.2. r— liminf inequalities. Next, we turn our attention to the proof of the sepa- 
rate F-hminf inequahties for energy and dissipation. Let us start with a statement 
concerning the energy densities. 

Lemma 3.2. Under assumptions (2.1e) and (2.6d), we have 

^ei ^ I ■ Ic ^'^^ ~^ I ■ Ih locally uniformly. (3.3) 
Moreover, we have 

\z\l < inf (hminf Wh'(^e) I z, ^ (3.4) 

L £— >0 ' J 



Proof. Let Kq ^ M'^^'*, fix (5 > and find the corresponding Cc\{S) > from 
condition (2.1e). As eKq C i?c^,(^)(0) for e sufficiently smaU we have that 

limsup sup IW^gi — I • |c| < <5 sup | • ^ < 5c 

e^O Ko Ko 

and local uniform convergence follows from 6 > being arbitrary. The same argu- 
ment apphes to W^. 

As for the F liminf inequality (3.4), let Zg ^ z and assume with no loss of 
generality that supgW^{zs) < oo. Hence, W^{zs) = W^(ze) and the inequality 
follows from the above proved uniform convergence. □ 

We are now in the position of proving the F-liminf estimate for the energy. It 
follows indeed from (3.3) and the lower-semicontinuity result of Lemma 4.2. 

Lemma 3.3 (F-liminf for the energy). For all {u,z) G Q we have 

Wo(u,z) < inf < liminf We («£, Ze) I (ucZs) [u, z) weakly in Q>. 

Proof. Let (?i£, Ze) — > (m, z) weakly in Q. Wc can assume with no loss of generality 
that sup^We{ue, Ze) < OO. Owing to the F-liminf inequality (3.4) and the lower 
semicontinuity Lemma 4.2 we readily conclude that 

/ b|^< liminf / z^) dx = liminf ^ [ Wh(I+eZe)dx. (3.5) 
Jn ^-^0 Jn ^ Jq 

Moreover, >Ve(we, z^) < oo implies ez^ £ K—I almost everywhere. In particular, 
eZe are bounded in L°°. The same holds for (J+ez^)^^ as 

(J+iTZe)"^ = C0f(/+£Ze)/det(/-t-£Ze) = Cof (I+eZs). 

We define the auxiliary tensors 

We := \ {^{l+eze)-^ -1 + eze) = e{I+ez,)-^z^, (3.6) 

so that (/+eze)~^ = I—eZg+eWs. By the first equality in (3.6) we have HeWellLo" < 
c, while the second gives 

llwelki - e||(/+eze)-iz2|lLi < ce\\z,\\l2 < ce 

where we have also used the boundedness in L^ of Ze from (3.1). Thus, by interpo- 
lation, We is bounded in L^ as well, so that Wg ^ weakly in L^. 
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Given Ae := (Fei,£— /)/e we want to show the weak convergence Ag — > \/u—z. 
Prom 

A, = -^ {{I+sVu,){I+sz,)-^ - I) (3.7) 
we find I+sA^ = {I+sVus){I+eZs)~^ and compute that 

Ae = ^{{I+eVUe){I-SZe+SWe) - l) = VWe - Ze + We - s{VUeZe-'VUe'We) . 

Hence, as we have that \/Ue—Ze Vu—z and We ^ weakly in L^, we have to 
show Ve := Vue{eze—ewe) — >■ weakly in as well. Indeed, the boundedness in 

of Ve follows from ||Vue||L2 < c (sec (3.1)) and the L°°-boundedness of eze and ewe- 
Moreover, since Ze and Wg are bounded in we have ||i'e||Li ^ ce and conclude 
Ve ^ weakly in L^. 

Eventually, owing to Lemma 3.2, we are in the position of exploiting the lower 
semicontinuity Lemma 4.2 in order to obtain that 

f |Vu-^|^ < liminf f W^i{Ae)dx 
Jq in 

= lim^inf ^ J Wei{{I+eyue){I+ezer^) dx. 

Finally, by recalling relation (2.3) and the already established (3.5) the assertion 
follows. □ 

Before moving to the F-liminf inequality for the dissipation functionals Ve, we 
prepare here a preliminary result on the functions De- 

Lemma 3.4 (F-convergence of De). De Dq in the sense of T -convergence. 

Proof. F-liminf inequality. Let {zej'Ze) — > (zyz) and assume with no loss of gen- 
erality that supg De{ze,Ze) < 00. In particular, we have that {I+eZe){I+eZe)~^ € 
SL{d). By defining 

Ce := J {{I+£%){I+SZe)~^ - I) =Ze - Ze + We - SZeZe + SZeWe 

where We is given in (3.6), we readily check that I+sCe € SL(rf) and (e ^ z — z. 
Let now t Pe{t) G C^(0, 1;K''^'') be such that Pe{0) = I, Pe{l) = I+eQ, and 



D{I,I+eQe)>{l-e) I R{PePr^)dt. 
Jq 



Such function Pe exists by the very definition of D. By possibly reparametrizing 
Pe and using assumption (2.10b) and bound (2.12) we can assume that 



(2.10b) . (2.12) 
c4Pe{t)Pe-\t)\ < R{Pe{t)Pe-\t))<2D{I,I+eCe) < CE. (3.8) 



Hence, 



\Pe{t)-I\ < ^ \PePe~'\ \Pe\ ds < CS \Pe\ds < CS (^1 + |P,-/| ds^ 

so that Pe ^ I uniformly by Gronwall Lemma. 
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By defining Pe{t) = 1+ {Pe{t)-I)/e one has that Pe{0) = I and Pe{l) = I+Cs- 
Moreover, as sPe = Pe and R is positively 1-homogeneous (2.10a), we have that 



1 A 
-D{I,I+eQ > {1-s) / R{P,p-^)At. 
^ Jo 



Owing now to bound (3.8), by possibly extracting not relabeled subsequences, we 
have that Pe ^ Q weakly-star in L°°(0, 1;M<*^'') and 

limmi De(zs,Ze) = liminf -D(I,I+e(^s) 

s—>0 e—>0 e 



> 



liminf/ R{P^p-'^)dt> [ R{Q)dt>R{Q) 
Jo Jo 



where we have exploited the lower semicontinuity tool of Lemma 4.2 and used 
Jensen's inequality with Q = Qdt. 

Finally, by integrating we have that 



z — z 



Q= f Qdt=\im f Pe di = hm 
Jo Jo 

so that we have checked 

limmf Ds{ze,%) > R{z—z). 

e— >0 

Recovery sequence. Given C, G ^"^^ we have that exp((^) e SL((i) and, by taking 
P{t) := exp(t<^) into the definition of D, we readily check that D{I, exp((^)) < -R(C)- 

Let now z G -I^(jgv given and define 

% = ^[e-x.-p{e{z-z)){I+ez) - l). 
As {I+£Zs){I+ez)~^ = exp(£(2"— ^;)), we have that 

limsup I?£(z, Je) = lim sup -£)(/, exp(e(z—2;))) < R{z—z) — Dq[z,z) 
so that [z, %) is a recovery sequence. □ 



Owing to Lemma 3.4, it suffices now to apply the lower semicontinuity result in 
Lemma 4.2 in order to establish the F-liminf inequality for the dissipation func- 
tionals. More precisely, we have following. 

Lemma 3.5 (F-liminf for the dissipation). 

Vo{z,z) < inil limini V^iz^^Zs) I 

{z„z,) {z,z) weakly in (L2(f2; M<'x'^))2|. (3.9) 
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3.3. Mutual recovery sequence. We now come to the construction of a mutual 
recovery sequence. Let us recall from [MRS08] that indeed two separate F-limsup 
inequalities for energy and dissipation generally do not suffice for passing to the 
limit in RIS. In particular, the construction of recovery sequences for energy and 
dissipation has to be mutually coordinated. 

Lemma 3.6 (Mutual recovery sequence). Let t e [0,T], {ue,Ze) — > {uq,zo) weakly 
in Q, and 

SUp£e{t,Ue,Ze) < 00. (3.10) 
s 

Moreover, let (uo,zo) := (mo,2o) + {u,S) with {u,z) e 0^(0; M*^) x C^(0; R^^^**). 
Then, there exist {ue,%) G Q such that 

(ug,%) — > {uo,zq) weakly in Q and 

lim sup (£e{t,Ue,Ze) - £e{t,Ue,Ze) +Vs{Ze,%)] 

< {£o{t,uo,%) - £o{t,UQ,zo) +Vq{zo,zo)^. (3-11) 

Proof. For the sake of clarity, we decompose this argument into subsequent steps. 
The general strategy of the proof is to choose (we,2e) and show convergence to 
(wo,2o), 

Y\Va.SW£i'Ds{Ze,%) < T>q{zq,Zo) = R{z), 

£->-0 

and 

lim sup (£e{t,Ue,Zs) - £e{t,Ue,Ze)] < £o{t,Uo,Zo) - £o{t,Uo,Zo). 

Note that in order to establish the latter we cannot argue on individual terms but 
rather aim at exploiting certain cancellations. This resembles the situation of the 
so-called quadratic trick (see, e.g., [MT05]) and crucially uses (2. Id) as well as the 
smoothness of (u,, z). In particular, note that within this proof the constant c may 

depend on u and z as well. 

Step 1: Choice of the mutual recovery sequence. By defining the functions tpe := 
id + eu and ip^ := id + eue and the set 

fie := {a; e ri | e^Tp{ez{x)){I+eZe{x)) e K^, 

the proof of the lemma follows by checking that the choices 

£ _ f ^{exj£>{ez){l+eze) - 1) on 9.^ 
\ Ze else, 

fulfill (3.11) and, in particular, [ue,Ze) — >■ {uq,zq) weakly in Q. The construction of 
Ue via a composition and oi'Ze via matrix exponential and multiplication is necessary 
in order to deal with the multiplicative nature of finite-strain elastoplasticity. 

Note that the construction of the mutual recovery sequence is compatible with 
the constraint det(/+£Vue) > considered in (2.1a). Indeed, by letting e be small 
enough we have that I+eVu is everywhere positive definite, hence det(/-|-£Vu) > 0. 
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In particular, as det(/+eVM£) > almost everywhere by (2.1a) and (3.10), we have 
that 

det(/+£Vu£) = det(VV'e(<^e)V</9£) = det{I+e'Vu{ips)) det{I+eVue) > 

almost everywhere as well. That is, I+eVue G GL+(d) almost everywhere. 

From the bound (3.10) we readily have that I+ez^ G SL{d) almost everywhere. 
Hence, upon noting that 



j^gg- = I exp{ez){I+ezs) on 
^ \ I+sZg else, 



we immediately check that {I+e%) ^ K C SL{d) almost everywhere and is bounded 
in L°° . Using the fact that tr z = wc have det exp{£z) = exp(£tr z) = 1 and hence 
exp(ei)(/+eze) e SL{d) almost everywhere. 

Next, note that the measure of the complement of can be controlled by means 
of a Chebyshev estimate. Indeed, relation (2.8) gives 

\n\ns\= ldx<cji \exp{ez){I+eze) - dx 

Now, one has that 

% - Ze = -(exp{ez){I+eze) - I) - Ze ^ -{exp{ez)-I){I+eze) on fl^, 
e e 

Zs — Zg = on f2\r2e, 

the convergence jfi \ fiej — > 0, and that % and z^ are bounded in L^. Hence, we 
readily check that 

%-Zs^z strongly in L^(f2; M'^^'^). (3.12) 
This implies that % ^ 'zq = zq+z weakly in L^, hence 

% + Zs^Zo + zo weakly in L^{n; M'^^''), . (3.13) 

Prom the energy bound (3.10) and the coercivity Lemma 3.1 we have that Ug 

is bounded in H"'^ and eu^ — > strongly in L^. Hence, one has that ||i^£— id||L2 = 
sll^ellL^ < c£ and, by the Lipschitz continuity of Vu, we conclude that 

\\Vu{'^^) - Vw||l2 < c||(/?£-id||L2 = C£||lte||L2 < C£. (3-14) 

Moreover, by computing 

S/Ue = i(VV'£(<^e)V<^e-/) = ^{{I+eyu){ipe)Vip,-l) 

we obtain that 

\\{Vu,-VUe) - Vti||L2 < ||Vw((/Je)-Vti||L2 + \\eVu{^p,)VUe\\i^2 
(3.14) 

< c£ + ce\\Vus\\i2 <ce (3.15) 
and this implies that Ug ^uq = uq+u weakly in H^. 
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The tensors Ae = (Fci,^-/)/e and Ae = {Fci,e-I)/e fulfill 

A, = i((J+£V«e)(/+e^e)-' - /), A, = l((7+eVue)(/+e2e)-' - /) 

and are hence both bounded in by (2.7). 

Fix now 6 and let Cc\{d) and Ch((5) be given by conditions (2.1e) and (2.6d), 
respectively. For all e > we define the sets 

:= {x e f2 I \eA,{x)\ + \eA,{x)\ < Cei(5)}, 

Z^^ := {a; e n I \ez,{x)\ + \ez,{x)\ < Ch{6)}, 

We refer to the latter as good sets as strains are there under control and we can 
replace the nonlinear densities Wei and Wh by their quadratic expansions via (2.1e) 
and (2.6d). In particular, on the good sets the quadratic character of the expansions 
will entail the control of the difference of the energy contributions by means of a 
suitable cancellation (quadratic trick). On the other hand, we term bad sets the 
corresponding complements f2 \ C/^ and fl\ where the quadratic expansions 
are a priori not available. Using some nontrivial cancellations, we will show that 
the difFcrcnce of the energy contributions on the bad sets is infinitesimal. Note 
preliminarily that the integrands on the bad sets blow up while the bad sets have 
small measure. Indeed, 

\n\u!\= [ ldx<-^ [ (\A,\ + \A,\fdx<-^ (3.16) 

P\Zl\= [ ldx<-^ [ i\z,\ + \%\fdx<-^ (3.17) 
Ja\zi Ci^(oj Jq Cjj(d) 

Step 2: Treatment of the dissipation term. As % = Ze on Q, \ fig one has that 

Vs{ze,%) = - [ D{I,exp{ez))dx <- [ D{I,exp{ez))dx. (3.18) 
^ Jne ^ Jn 

In the construction of the recovery sequence in the proof of Lemma 3.4 we have 
proved that 

lim sup -D{I, exp{sz)) < R{z). (3.19) 

Eventually, by taking the lim sup in relation (3.18) and using (3.19) we have proved 
that 

lim sup T>^{zs, %) = lim sup - f D{I, exp{ez)) dx 

e^o £-5-0 e Jn 

< [ R{z)dx = Vo{zo,zo). (3.20) 
Jn 

Step 3: Limsup for the differences of the elastic energy terms. Let us start by 
rewriting the tensors A^ as 

As = WUs - Zs + We - sVUeZe + eVUsWe 
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where We is given by (3.6). On the other hand, as regards the tensors we have 
that 

Ae = ^ {{I + eVue){I-sZe-\-eWe) exp(-ez) - I) 

= {VUe-Ze+We-eVUeZe + sVUeWe) exp(-e5) 

+ -(exp{—ez)—l) on fig 

Ae = l{i^ + £'^Ue){I-eZe+eWe) - l) 

= Vus—Ze+Ws—sVusZc+eS/UsWc on Oyn^. 
Hence, one can compute that 

A^-A^ = {Vus-Vue){I-ez^+swe) + ^{cxp{-eS)-l) 

+ (Vm£— eVM£Z£+£VM£U>£)(cxp(— ez)— /) on fig 
Ae- A, = {yue-Vue){I-eze+ewe) on Q\ Q^- 

In particular, owing to convergence (3.15) and the L°° bounds for ezg and eWe (see 
the discussion after (3.6)) we have that {'Vus—'Vue){I—eZe+eWs) converges to Vu 
strongly in . Thus, by recalling that — > weakly in it is a standard matter 

to check that 

Ae + Ae^ (Vwo-zo) + (Vmo-^o) weakly in L2(f]; R''^''), (3.21) 
Ae-Ae^Vu-z strongly in L^{n; W^""^). (3.22) 

On the good set we will use the assumption (2.1e) in order to have that 

W^Me) W^,{Ae) < - \Ae\l + 6{\Ae\l + 
= ^(Ae-A,y.C{A+Ae) + 5{\A,\l + \A\l). (3.23) 

Let us now argue on the bad set O \ f7* by defining 

Gi^e := {I+e\7u,){I+e\/ue)-\ Gz.e := (/+eze)(/+£^e)"^ 

The energy bound (3.10), together with assumption (2.1a), implies that I+eVUg 
is invertible almost everywhere. Note that and G2,e are chosen in such a way 
that Fei,e = Gi^eFei,sG2,£- We readily compute that 

Gl,e -1 = VV'e(V'£)V</Je(/+£VUe)-^ - / = ^^MVe) -1 = eVu(v>e) 

SO that ||Gi,e-/||L~(n\[/|;R''xd) = s\\Vu{'fis)\\i,,^((i\us.]gtdxd) < cs. Moreover, one 
has that 



exp(-e^) on {n \ U^) n Vl^, 
on f]\ (C/^* UQe). 



Hence, \\G2,e — — as well. Next, estimate (4.1) and bound (3.10) 

allow us to control the elastic part of the energy on the bad set fi \ f/^ (where Vue 
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and Ze are not under control) by cancellation. For this we employ the multiplicative 
estimate (2. Id) provided in (4.1): 



£■ 
(4.1) 



\ I (W;i(Gl,eFel,eG'2,.) - W,x{F,x,e)) dx 
'< % [ (W^el(^^el,e) + Cg) (|Gl,e-/| + |G2,e-/|) dx 

< C7 (^^ W^cl(FeLe) dx + U',]^ (||Gi,e-J||L~ + ||G2,e-/||L~) 

(3.10)&(3.16) / 1 \ 

< c(l+g^). (3.24) 

Thus, we have controlled the difference of the energy contributions in the bad set 
n \ where the gradients are big. 

Finally, by using convergences (3.21)-(3.22), equation (3.23) on the good set , 
relation (3.24) on the bad set Q \ U^, and the boundedness of and A^, we 
conclude that 



limsupf/ W^,{A,)dx- [ W^,{A,)dx^ 
e^O \Jn Jq j 

(3.23) / 1 /■ 

< limsup - / (Ae-Ae):C(Ae+A,)dx + C(5 

£-)-0 \^ Jul 

+ / iWlU,) - Wl.^A,)) dx] 

(3-24) ( I r ^ ^ ( 

< limsup - / {A^-A^-£.{A^+A^ da; + c5 + c 1-^ 

e^O y 2 \ 

<\l (Vu - 5) :C(V(uo+wo) - (2o+-2o)) da; + d> 

-sym csymi2 / i^y., sym .^symi2 



V«;;^"'"-2^^'"|^dx^ / |V<'"-z^''"^|^dx + c^ (3.25) 
Jn 

where we have made use of relation (2.3). 

Step 4-' Upper bound on the hardening energy term. Let us now turn our attention 
to the hardening part of the energy. On the good set we have that 

= ^{%-z,yM{%+z,) + S{\%g+\z,g). (3.26) 
As regards the bad set i^\Z^ one has that 

Wl{z,)-W^{ze) = I ^W?h(exp(£5)(/+£Ze)) - ^Wi,{I+ez,) on {n \ Zl) n 

" ^ ^ \ onO\(ZfuOe)- 
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Hence, by exploiting the local Lipschitz continuity of Wh we have that 

[ {Wl{%) - W^{z,)) drr < 4 / I exp(£5) - I\ \I+eZe\ dx 
Jn\zi £ Jn\zl 

Eventually, owing to convergences (3.12)-(3.13) we compute that 
limsup f/" Wl{%)dx~ [ Wl{ze)dx] 

(3.26) ( f 

< limsup I j -{ze-Ze)-M(^s+Ze)dx + c5 



>0 



{W^{%) - Waze)) dx 



ce 



(3.27) / r 1 

< limsup / -{%-z^)-M{%+Zs)dx + cS + ^— , 

= [ ^z-M{zo+zo)dx + cS= [ l^oliida;- / \zo\idx + c6. (3.28) 
Jn ^ Jn Jn 

Step 5: Conclusion of the proof. By collecting relations (3.25) and (3.28), and 
recalling that {£{t),Ue — Us) {^(t), uq — uo) we have proved that 

]imsup (^Ss{t, Us, %)-£s{t, Us, Zsj^ < (^£o{t,uo,zo)-£o{t,uo,zo)^ +cS. 
Finally, the assertion (3.11) follows by taking 5 — > and employing (3.20). □ 



3.4. Proof of Theorem 2.2. Owing to the the above-obtained F-liminf and 
mutual-recovery-scqucnce results, the proof of Theorem 2.2 now follows along the 

lines of the general theory of [MRS08] . Wc limit ourselves to sketch the main points 
of the argument and refer the reader to [MRS08] for the details. 

Let {ug,Zs) be a sequence of finite-plasticity solutions. The coercivity of the en- 
ergy (3.1) entails an a priori bound on {u^, z^). In particular, we have the following. 

Corollary 3.7 (A priori bound). There exists c > such that all finite-plasticity 
solutions {Ue,Ze) fulfill 

Vi G [0,T] : ||Ue(t)||Hi + |ke(t)||L2 + ||£^e(t) ||loc + DisSp, (^e! [0,*]) < C. (3.29) 
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Proof. We exploit the energy balance (2.16) and the bound (3.1) in order to get 
that, for all t G [0,T], 



\\Vu,{t)\\l. + \\z,{t)\\l. + \\ez,{t)\\loo+BissT,^{ze; [0,t]) 

(3.1) 

< c{l+We{u,{t),Z,{t)))+msST:>AZe;[0,t\) 

< c(l + £e{t, U,{t), Ze{t)) + {l{t), Ue{t)) + DisSp, {Ze', [0, t])) 

< C [l + £,%U%zl) + \\i{t)U-A\uS)W + ^ ll^llH-ill^^ellHi ds 



so that the assertion follows by Gronwall Lemma. □ 



Owing to the a priori bound (3.29), we may now exploit the generalized version 
of Helly's Selection Principle in [MRS08, Thm. A.l] (consider also the comments 
thereafter) and deduce that, at least for some nonrelabeled subsequence, and all 
s, t e [0, T] with s<t, 



So{t) limDissi,.(ze;[0,t]), 

£—7-0 

Ze(t) — >■ ZQ{t) weakly in Z, 
Diss2?o(2o; [s,t]) < 6o{t) - 6o{s), 



Moreover, by letting t £ [0,T] be fixed we may extract a further subsequence 
(still not relabeled, possibly depending on t) such that Ue{t) — >■ weakly in U. 
We now check that indeed (u,, zo{t)) S So{t). To this aim, by density it suffices to 
consider competitors (^05^0) = {u^,, Zo{t)) + {u, z) with {u,z) smooth and compactly 
supported. By applying Lemma 3.6 we find a mutual recovery sequence (ug, %) such 
that 



£o{t,uo,zo) - £oit,u^,ZQ{t)) +'Do{zo{t),za) 

> limsup {£,{t,Ue,Ze) - £s{t,U,{t),Z,{t)) + 'De{Zs{t),Ze)) > (3.30) 

£->-0 



where the last inequality follows from the stability (2.15) of {u^{t), z^{t)). Hence, we 
have proved that (w*,2:o(i)) € So{t). Note that, given zo{t) G Z, as the functional 
u G U £o{t,u,ZQ{t)) is uniformly convex there exists a unique Uo{t) G U such 
that {uo{t), zo{t)) G So{t). From the fact that (u*,zo(i)) G 'Sa{t) we conclude that 
= Uo{t). In particular Ue{t) uo{t) weakly in U for all t G [0, T] and the whole 
sequence converges. 
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Let now be given a partition {0 = to < ti < ■ ■ ■ < In = t} . By passing to the 
liminf in the energy balance (2.16) and using Lemmas 3.3 and 3.5 we get that 



N 



So {t, Uo {t) ,Zo{t))+^ Vo {Zo {U ),Zo{ti-i)) 

< lim^inf (^£e{t,Ue{t),Ze{t)) + Y,V,{z,{ti),ZeiU-i))^ 

< liminf f£:,(0,uO,zO)- f {£,u,)d.s) = So{0,u°, z^) - f {i,uo)ds 







where for the last cquaHty we have used (2.14) and the convergence of u^. Hence, 
the upper energy estimate follows by taking the sup among all partitions of the 
interval [0,t]. The lower energy estimate can classically recovered from stability as 
in [Mie05, Prop. 2.7]. This proves that (uq, Zq) is a linearized-plasticity solution. In 
particular, as linearized-plasticity solutions are unique, the whole sequence (w^, Zg) 
converges and no extraction of subsequences is actually needed. 

Along the lines of the proof of Theorem 2.2 (see also [MRS08, Thm. 3.1]) we 
also obtain the following convergences. 

Corollary 3.8 (Improved convergences). Under the assumptions of Theorem 2.2 

we have that, for all t G [0,T], 

[ {W:,{A,)+Wl{z,))dx^ [ {\Vuo-zo\l + \zog)dx, (3.31) 
Jn Jn 

Dissv,{ze;[0,t]) ^ R{z)ds. (3.32) 

^0 

In particular, owing to the energy convergence (3.31) we are in the position of de- 
ducing some strong convergence of finite-plasticity solutions to linearized-plasticity 
solutions. 

Corollary 3.9 (Strong convergence). Under the assumptions of Theorem 2.2 for 
all t G [0,T] we have that {ue{t),Ze{t)) -)■ {uo{t),zo{t)) strongly in W^'P(0;M'') x 
LP(f2;M''x'') for all p & [1,2). 



Proof. Let v denote the Young measure generated by the sequence (^e,Ze) and 
define the measure (^^Z) := j/(A''®Rf^^f , Z) for all Borel sets (A*,Z) C 
K^y^^ X R''^''. Note that i/"*"" is indeed the Young measure generated by {Af^, z^). 
By using the lower semicontinuity Lemma 4.2 and the energy convergence (3.31) 
we deduce that 

{\A^y^\l + \z\l)dvT\A^^^,z) \ dx 



< 



{\A\l + \z\l)d,y,{A,z)) dx 
liminf / {Wl^{A,) -h Wi{z,)) dx ^ I (|Vuo-^o|^ + \zo\V) (3-33) 
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Recall from (3.7) that 

where the remainder term e{'VugZg~Vu^'w^Y^'^ converges strongly to in L^' for 
all p e [1, 2). Hence, the baryccnter of j/^^y™ is clearly (Vmq*'™— Zq^'", zq). 

We readily check that the measure i/'^^™ is concentrated in its barycenter. Indeed, 
if this was not the case, by uniform convexity we would have that 

/ {\Vur-zr\l + \zog)dx 

< f i , (|A«^-|^ + |z|^)d^f-(A«''-,^)] dx 

contradicting relation (3.33). Here we have used positive definiteness from (2.5) 
and (2.9). As v^'^'^ is concentrated, we exploit [AGS08, Thm. 5.4.4.iii, p. 127] and 
deduce that 

/ f{x,Af^{x),z,{x))dx^ I (I f{x,A^y^,z)dvf^{A^y^,z) \ dx 

Jn Jq. yiR^yl^xM<*x<* j 

along with the choice 

f{x,A^y'^,z) := \{Wu''r{x)-z''nx),z''^'^{x)) - {A^y^-z^y^,z)\r 

Hence, we have that {Af^.z,) (Vm^'''"-^^^'", 0o) strongly in Lf(0;Kfy><^) x 
U>{n- W^^^d.) for all p e [1, 2). In particular, 

V«7'" = A7'° + 0|'"" + £(Vue^,-V«e«'e)''""^ Vmo'"" strongly in 
for all p G [1,2) and the assertion follows by Korn's inequality. □ 

3.5. Sketch of the proof of Theorem 2.3. The argument for Theorem 2.2 can 
be adapted to prove Theorem 2.3 as well. The only notable difference is that one 
has to cope with the fact that the piecewise constant interpolants {lle,Ze) of the ap- 
proximate incremental minimizers need not be stable but rather just approximately 
stable. More precisely, from (2.17) and the triangle inequality we have that 

V(w,t) G Q : £S^u,z) - £S^uS),zS)) + > -Tea^. 

By coordinating to the sequence {u^{t),'Zg(t)) a mutual recovery sequence (ug,?^) 
via Lemma 3.6 (with {us{t),Zs{t)) instead of {ue{t), Zs{t))) the lower bound (3.30) 
still follows as Tgag — > 0. Hence, the stability of the limit can be recovered. Fi- 
nally, improved and strong convergences in the spirit of Corollaries 3.8-3.9 can be 
established as well. 

Acknowledgement. We are gratefully indebted to the Referee for her/his careful 
reading of the manuscript. 
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4. Appendix 

4.1. Estimate on left and right multiplication. In the proof of Theorems 2.2- 
2.3 we have made use of the following estimate combining left and right multipli- 
cation. 

Lemma 4.1. Assume (2.1a) and (2. Id). Then, 

3c7, C8,7 > VGi, G2 e B^{I) VF e GL+{d) : 

|Wel(GlFG2)-Wel(i^)| <C7(W^(F)+C8)(|Gl-/| + |G2-/|). (4.1) 

Proof. Following [Bal02, Lemma 2.5], we find positive constants cq, cq, 7 such that, 
for all G e B-f{I) and all F e GL+(d), one has that 

WeliGF) < CoWel(F) + CO, Wel(FG) < CoWel(F) + CO, (4.2) 

\dFWiGF)F^\<coW,,iF)+co, (4.3) 

|i^TaFW^(i^G)| < coVFci(F) + CO. (4.4) 

For s G [0, 1], let now Hj(s) := (l-s)/+sGj for j = 1, 2, and note that Hj e -B7(/). 
As the derivative Hj = Gj — I is constant we can compute that 

W,i{GiFG2)-W,i{F)= —W,i{H,{s)FH2{s))ds 

= [ dFW^i{HrFH2):{H[FH2+HiFH'2)ds 
Jo 

= [ dFW^x{HrFH2){FH2y ds : H[+ [ {H^Fy dFW^x{HrFH2)ds : H'^. 
Jo Jo 

We control the above right-hand side as 

^ dFW^i{HtFH2){FH2)^ ds : H[ (^j^ {coW^i{FH2)+co) d^\Gr-I\ 



^<'^ (c^Wei(F)+coCo+co)|Gi-7|, 



1 

[HiF)^ dFW^i{HiFH2) ds : H'^ 





(4.4) 
< 



^ (coWel(i?lF)+Co) d^ \G2-I\ 



^<'^ (c^W^el(i^)+C0C0+C0)|G2-/|, 



whence the assertion foUows. □ 

4.2. Lower semicontinuity tooL In Section 3 the following lower-semicontinuity 
lemma is used. 

Lemma 4.2 (Lower-semicontinuity). Let fo, fe : R" — >■ [0,c«] be lower semicon- 
tinuous, 

Vuo e M" : fo{vo) < inf { li^nf /^(ue) | -)■ t^o}, 

and We — >■ Wo weakly in L^(f2;]R"). Denoting by v the Young measure generated by 
We we have that 

/ ( / fo{w)di'a:{w)] da;<liminf / fs{ws)dx. 
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In particular, if fo is convex we have 

/ /o(wo) da; < liminf / fe(we)dx. 
Jn Jn 

This lemma is in the same spirit of the results by Balder [Bal84a, Thm. 1] 
and lOFFE [Iof77] and can be proved via augmenting the variables by including 
the parameter e. The reader is referred to [Ste08, Thm 4.3, Cor. 4.4] or [MRS09, 
Lemma 3.1] for a proof in the case d = 1. In case of local uniform convergence, a 
proof can be found in [Li96] . 
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